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I .  INTRODUCTION 


Over  the  past  several  years  a  number  of  schemes  have  been  devised 

for  performing  useful  optical  signal  processing  functions  in  an  integrated 

(12) 

optics  format.  *  One  of  the  most  attractive  of  these  is  an  acoustooptic 
spectrum  analyzer^  now  being  developed  under  Air  Force  sponsorship.  In  the 
device  a  waveguided  beam  is  Bragg-diffracted  by  a  surface  acoustic  wave.  The 
diffracted  beam  is  focused  by  a  lens  on  a  detector  array,  where  there  exists 
a  one-to-one  correspondence  between  the  element  in  the  array  that  is  illumi¬ 
nated  and  the  oscillation  frequency  of  the  surface  acoustic  wave. 

One  problem  yet  impeding  the  development  of  a  versatile  spectrum 
analyzer  is  optical  waveguide  scattering.  Owing  to  the  propagation  of  the 
processing  beam  parallel  to  and  just  beneath  the  waveguide  surface,  scattering 
centers  associated  with  surface  roughness  and  defects  in  the  top  one  micron  of 
material  have  a  devastating  effect  on  the  information  content  of  a  guided  wave. 
This  is  true  even  for  surfaces  and  materials  of  such  high  quality  that  a  beam 
traveling  normal  to  the  surface  is  perturbed  by  an  almost  immeasurable  amount. 

An  example  showing  the  deleterious  effects  of  optical  waveguide 
scattering  on  spectrum  analyzer  performance  is  that  in  which  two  closely  spaced 
frequencies  are  simultaneously  present  on  the  surface  acoustic  wave,  one  of 
them  in  much  greater  strength  than  the  other.  In  this  event,  scattered  light 
from  the  stronger  dif fracted-beam  component  can  arrive  at  the  same  detector 
element  as  the  unscattered  weaker  diffracted  beam-component.  The  weaker  signal 
will  be  partially  or  totally  masked.^  The  effect  is  one  that  limits  the  dynamic 
range  of  the  system.  To  avoid  the  problem,  small-angle  in-plane  waveguide 
scattering  must  be  reduced  in  intensity  to  the  point  where  other  sources  of 
spurious  signal,  such  as  those  resulting  from  lens  aberrations  and  positioning 
errors,  dominate.  As  integrated  optics  component  technology  improves,  increas¬ 
ing  emphasis  will  be  placed  on  the  need  for  reducing  waveguide  scattering  to 
the  lowest  possible  levels.  This  has  been  the  goal  of  the  present  program. 

In  this  report  we  describe  experimental  and  theoretical  approaches 
used  toward  reaching  this  goal,  as  well  as  insight  that  these  approaches  have 
provided  regarding  the  physical  nature  of  the  scattering  mechanism.  Waveguides 
formed  by  Ti  diffusion  into  LiNbO^  substrates  have  been  emphasized,  since  these 
are  the  current  waveguides  of  choice  for  the  spectrum-analyzer  application. 


However,  many  of  the  methods  and  procedures  are  applicable  to  other  types  of 
waveguides,  including  thin-film  waveguides  sputtered  on  S102-S1  substrates. 

The  achievements  of  the  program  reported  here  include  the 

following: 

•  A  theoretical  analysis  of  both  in-plane  and  out-of-plane 

(3) 

scattering  that  extends  the  work  of  Marcuse  and  Boyd 

(4) 

and  Anderson  ; 

•  A  modulation  technique  for  studying  weak  out-of-plane 
scattering  in  the  presence  of  scattering  from  non-waveguide 
sources; 

•  A  50  dB  dynamic  range  system  for  studying  in-plane  scatter¬ 
ing; 

•  Controlled  polishing  techniques  for  removing  scattering 
material  at  or  near  the  surface  of  an  optical  waveguide; 

•  Demonstration  of  the  utility  of  Nomarski  microscopy  for 
examining  the  surface  of  optical  waveguides  and  substrates; 

•  Acquisition  of  in-plane  scattering  data  at  several  visible 
wavelengths ; 

•  Development  of  novel  diagnostic  methods  for  obtaining  infor¬ 
mation  about  sources  of  scattering,  including  methods  based 
on  electrooptic,  photorefractive,  and  superstrate  refractive- 
index  effects. 

These  achievements  have  enabled  us  to  draw  the  following  conclusions  regarding 
the  nature  of  scattering  in  Ti-diffused  LiNbO^  waveguides: 

•  Bulk  scattering,  at  least  some  of  which  is  diffusion-induced, 
is  far  more  important  than  surface  scattering. 

•  The  bulk  scatterers  are  anisotropic  in  refractive  index  and 
possibly  in  shape. 

•  The  scatterers  are  associated  with  regions  that  have  different 
electrooptic  properties  from  those  of  the  crystal  as  a  whole. 

•  Scattering  increases  in  proportion  to  the  amount  of  Ti  used 
to  make  the  waveguide,  but  is  not,  within  limits,  sensitive 
to  the  time  and  temperature  of  the  diffusion  treatment. 

•  Scattering  is  similar  for  x-cut  and  y-cut  crystals. 


2 


The  recipe  that  we  employ  for  our  best  waveguides  requires  the  use 

O  0 

of  thin  150  A-300  A  Ti  layers  diffused  for  3  h  at  950°C  in  flowing  O2,  followed 
by  a  rapid  quench  to  600°C  and  a  slow  cool  to  room  temperature.  For  waveguides 
formed  using  thicker  Ti  films,  some  improvement  in  quality  is  always  obtained 
by  lightly  polishing  the  finished  waveguide,  even  when  all  Ti  has  been  dif¬ 
fused  into  the  surface. 

The  above  recipe  applies  for  waveguides  used  with  beam  propagation 
perpendicular  to  the  c-axis.  For  the  parallel  configuration,  however,  our 
best  low  scattering  performance  was  obtained  using  a  thick  Ti  film  (about 

O 

720  A)  and  a  long  post-fabrication  polish — long  enough,  in  fact,  to  remove 
most  of  the  initially  formed  waveguide  layer. 
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II.  THEORY  OF  OPTICAL  WAVEGUIDE  SCATTERING 


Most  of  the  data  presented  in  this  report  show  the  angular  distribu¬ 
tion  of  light  scattered  from  an  initially  collimated  waveguided  beam.  Since 
individual  scattering  imperfections  are  not  observed,  only  indirect  informa¬ 
tion  about  their  nature  is  available.  The  accuracy  and  extent  of  this  infor¬ 
mation  depends  on  the  realism  of  physical  models  used  to  describe  scatterers 
and  on  the  completeness  of  the  theory  employed  to  correlate  scattering  data 
with  these  models.  Accordingly,  considerable  effort  has  been  expended  on 
theoretical  aspects  of  the  waveguide  scattering  problem,  the  results  of  which 
are  described  in  this  section. 

Two  approaches  to  the  problem  are  characterized  as  the  Green' s-f unc¬ 
tion  approach  and  the  coupled-mode  approach.  The  former  was  implemented,  for 
example,  by  Giallorenzi,  ^  and  the  latter,  by  Marcuse^  .  When  carried  out 
with  complete  rigor  both  should  give  equivalent  descriptions  of  the  waveguide 
scattering  problem.  It  is  our  belief  that  the  Green's  function  approach  is 
more  suitable  for  obtaining  the  scattered  energy  distribution  from  an  indivi¬ 
dual  imperfection  of  specified  shape.  However,  if  this  is  not  required,  the 
coupled-mode  approach  is  more  suitable  since  it  allows  for  convenient  statis¬ 
tical  averaging  over  an  ensemble  of  scattering  centers,  each  of  which  may  be 
slightly  different. 

(3) 

We  concentrated  on  the  coupled-mode  formalism  for  two  reasons: 
first,  it  is  mathematically  simpler;  second,  the  results  are  directly  appli¬ 
cable  to  experimental  scattered  energy  distributions. 

THE  SCALAR  WAVE  EQUATION 


A  classic  exposition  of  the  coupled-wave  formalism  as  it  applied  to 
optical  waveguides  is  found  in  the  book  Light  Transmission  Optics,  by 

(3) 

Marcuse.  He  applies  the  formalism  to  the  particular  problem  of  scattering 

by  surface-roughness  fluctuations  in  symmetric-slab  waveguides.  The  surface 
roughness  is  assumed  to  vary  only  in  the  direction  of  beam  propagation.  For 
these  reasons,  the  theory  is  not  directly  applicable  to  LiNbO^  and  sputtered 
thin-film  waveguides,  which  are  asymmetric  rather  than  symmetric,  and  which 
have  volume  as  well  as  surface  scattering  centers  that  are  three 
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dimensional.  Nevertheless,  modifications  required  to  generalize  Marcuse's  re¬ 
sults  are  relatively  straightforward.  In  point  of  fact,  Walter  and  Houghton 
extend  Marcuse’s  theory  to  the  case  of  the  asymmetric  slab  waveguide, ^  and 

Boyd  and  Anderson,  using  heuristic  arguments,  extended  it  to  include  two  di- 

(4) 

mensional  scatterers.  Our  interest  is  in  (1)  providing  a  mathematical  basis 

(4) 

for  the  conclusions  of  Boyd  and  Anderson  and  (2)  extending  the  work  of  all 
authors  mentioned  to  include  the  effects  of  volume  scattering  centers  as  well 
as  surface  roughness. 

The  geometry  for  the  calculation  is  shown  in  Fig.  1:  x  is  the  direc¬ 
tion  of  propagation,  z  is  the  axis  normal  to  the  waveguide,  and  0  and  4>  are  the 
polar  and  azimuthal  angles  defined  in  the  usual  way.  For  a  given  power  in  the 
(assumed)  single  waveguide  mode,  we  calculate  the  out-of-plane  power  radiated 
into  the  solid  angle  d4>d8sin0,  as  well  as  the  in-plane  power  scattered  into  d$. 
The  electric  field  amplitude  E  is  governed  by  the  scalar  wave  equation, 


E  +  E  +  E  +  k2 

xx  yy  zz  o 


N2(x,y,z)£ 


0 


(1) 


2  2 

N  (x,y, z)  =  n  (z)  +  n(x,y,z) 

where  kQ  is  the  magnitude  of  the  free-space  wave  vector,  N  is  the  total  refrac¬ 
tive  index,  n  is  the  refractive  index  of  the  unperturbed  waveguide,  and  n  is  the 
refractive-index  perturbation.  We  will  consider  an  unperturbed  refractive- 
index  profile  like  that  shown  in  Fig.  2,  having  substrate  index  n2»  wave¬ 
guide  index  n^,  and  superstrate  index  nQ.  We  assume  that  the  substrate  index 
n^  is  very  close  to  the  guide  index  n^.  This  is  a  good  approximation  in 
LiNbO^  and  it  greatly  simplifies  the  mathematics. 

For  bulk  scatterers  having  an  average  rms  index  change  |6n|,  we  model 
the  refractive  index  variation  using  the  equations 


n(x,y,z)  =  n(z)g(x,y)  (2) 

n(z)  =  2n2<5n  z<0 

n(z)  =*  0  z >0 

That  is,  the  guide  index  profile  is  like  that  indicated  by  the  dashed  line  in 
Fig.  2,  where  the  position  of  the  dashed  line  relative  to  the  unperturbed 


Fig.  1.  Geometry  for  waveguide  scattering 
calculations . 


C 


n2(z) 


Fig.  2.  Refractive-index  profiles.  Solid  line:  Unperturbed 
profile.  Dashed  line:  Profile  in  the  vicinity  of  a 
volume  scattering  center. 
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solid-line  profile  varies  over  the  waveguide  surface  according  to  the  value 
of  g(x,y).  By  neglecting  the  z-dependence  of  6n,  we  assume  that  the  scatter¬ 
ing  centers  are  uniform  over  depths  at  least  as  large  as  the  waveguide  mode 
confinement  depth. 

(3) 

For  scattering  from  surface  roughness,  we  follow  Marcuse  in 

writing 

2  2 

n(x,y,z)  =  (nx  -  nQ)f(z),  (3) 

f(z)  =  +  1  in  (0<z<oog(x,y>) 
f(z)  =  -  1  in  (oog(x,y)<z<0) 
f(z)  =  0  elsewhere 

where  is  the  rms  surface  roughness. 

For  both  surface  and  volume  scattering  centers,  the  random  function 
g(x,y)  describes  the  spatial  distribution  of  scattering  centers  over  the  wave¬ 
guide  surface.  At  any  point  (x,y)  the  value  g(x,y)  indicates  either  the  index 
perturbation,  relative  to  6n,  or  the  surface  height  perturbation,  relative  to 
oq.  Thus,  g(xo,yo>  =  0  means  that  there  is  no  perturbations  at  point  (xQ,yo), 
while  g(x  ,y  )  =  1  means  that  the  perturbation  has  its  rms  average  value  at 

(Vyo>* 


SOLVING  THE  SCALAR  WAVE  EQUATION 

A  solution  to  the  scalar  wave  equation,  Eq.  (1),  is  obtained  by 
expanding  the  electric  field  in  terms  of  the  discrete  guided  modes  and  the  con¬ 
tinuum  of  radiation  modes.  The  single  discrete  mode  of  the  waveguide  we  are 

considering  has  a  spatial  variation  of  the  electric  field  that  may  be  written 
i B  *r 

EQ(z)e  — o  — ,  where  Eq(z)  is  the  modal  amplitude  distribution  normal  to  the 
waveguide  plane,  6^  *  0Q(c osiji,  sin<j>,  0)  is  the  wavevector,  and  £  =  (x,y,0)  is 
the  position. 

Each  element  of  the  continuum  of  radiation  modes  has  a  field  distri¬ 
bution  that  may  be  similarly  written  E(p,z)e*—  — ,  where  E(p,z)  is  the  modal 
amplitude  variation  of  the  radiation  mode  characterized  by  the  continuous 
parameter  p,  and  _B  is  the  wavevector  analogous  to  the  guided  mode  wavevector 
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that  describes  the  variation  of  phase  of  the  radiation  mode  in  the  waveguide 

22  2  1/2 

plane.  Each  radiation  mode  is  described  by  a  real  value  of  p  =  (k^  ) 

in  accord  with  Maxwell's  equations  and  appropriate  boundary  conditions.  De¬ 
rivation  of  the  forms  of  Eq(z)  and  E(p,z)  are  given  in  Ref.  (3)  for  symmetric 
waveguides  and  in  Ref.  (5)  for  asymmetric  waveguides.  Also  presented  in  these 
references  are  the  orthogonality  relationships  satisfied  by  the  various  dis¬ 
crete  and  radiation  modes,  which  we  will  later  have  occasion  to  employ.  In 
particular,  we  note 


Prad<5(p-p')  =  6 / 2 u3U o  /  E(p  ,z)E*(p'  ,z)dz 


(4) 


0  =  6o/2o3po  f  Eq(z)E  (p,z)dz 


P  =  SQ/ 2uuo  /  |Eq(z)|  dz 


In  these  expressions,  P  ^  is  the  power  in  the  radiation  field  contained 
within  an  infinitesimal  range  of  the  parameter  p,  and  P  is  the  power  in  the 
waveguide  mode.  Both  powers  are  given  for  beams  of  unit  width. 

In  the  case  of  weak  scattering  of  an  initially  guided  beam  propagat¬ 
ing  in  the  x-direction,  we  may  define  scattering  functions  cq((J>,x)  and 
q(p,<t>,x)  that  describe  the  amplitude  of  the  field  scattered  at  an  angle  4>  into 
the  guided  and  radiation  modes,  respectively.  In  terms  of  these  functions,  the 
total  field  may  be  written 

Ztr 

E(x,y,z)  =  Eo(z)ei^ox  +  j  d<|>co((l>»x)E0(z)ei^o^  —  (5) 

o 

2<|y  00 

+  /  d4>  /  dp  q(p,<(),x)E(p,z)ei-^^  —  , 

o  ° 


^  C<f>)  =  Bo(cos<t>,sini{>,0)  , 

1(41)  =  0^2  -  P2)1</2  (cos<t>,sin$,0)  , 

r  *  (x,y,0) 


9 


The  quantities  cq  and  q  are  determined  by  inserting  Eq,  (5)  into 
Eq.  (1)  and  making  use  of  the  orthogonality  conditions  given  in  Eq.  (4).  If, 
for  example,  both  sides  of  the  equation  obtained  by  inserting  Eq.  (5)  into 
Eq.  (1)  are  multiplied  by  EQ(z)e  •Eovv  and  the  resulting  expressions  are 
integrated  over  z,y,  and  <j>,  we  obtain 


x 

0  «  (2ir/6o|cos<j> 1  |  )  /  dx'rc”(<fr\x')  +  2iBocos*>’ 

o 


C^\x')]  +  k*  no  G( <J> '  ,x) 


G(<j>'  ,x) 


rx  ,  ,  (  ,  ,  ,  ,  iB  x  -  iB  ($') *r 

j  dx’  J  dy  g(x  ,y  )e  o  -o  — 

o  beam 

width  W 


(6) 


The  term  n  is  defined  as  the  part  of  the  integral  (B  /2wp  P)  /  dz  n(x,y,z) 
|E  (z)|  that  is  proportional  to  g(x,y).  Thus,  depends  on  the  model  for 

scattering  center  employed. 

A  result  similar  to  Eq.  (6)  is  found  for  the  radiation-mode  ampli¬ 
tude  q  when  we  insert  Eq.  (5)  into  Eq.  (1),  premultiply  both  sides  by 
E(p',z)e  ^  ,  and  integrate  over  z,y,<j>  and  p.  The  result  can  be  written 

using  Eq.  (6)  by  making  the  substitutions 

co(<t>',x')  -*  q(p'  ,<t>'  ,x)  ,  (7) 

B^*’)  £(♦’) 

and  by  redefining  nQ  to  mean  the  part  of  the  expression  (8/2wpoPratj) 

.oo  * 

J  dzE  (p',z)E  (z)n(x,y,z)  proportional  to  g(x,y). 

— oo 

Equation  (6)  for  cq  or  its  counterpart  for  q  are  most  simply  treated 

by  assuming  that  the  functions  c  and  q  change  slowly  enough  that  their  second 

0  (3) 

derivatives  may  be  neglected.  Then,  following  Marcuse  ,  we  consider  for¬ 
ward^)  and  backward (-)  traveling  waves  satisfying  the  boundary  conditions 

c£+)(*\0)  =  0  =  q(+)(p\$’,0)  ,  (8) 

c(_)(<fr\L)  =  0  =  q(-)(p\*\L) 
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where  L  is  the  length  of  the  scattering  medium,  measured  in  the  direction  of 
propagation.  The  slowly-varying  solutions  to  the  coupled-wave  equations  that 
satisfy  these  boundary  conditions  are  simply 

c£+)(*,L)  =  c£-)(*,0)  =  -k^(no)inG(4>,L)/4TTi  (9) 

q^(p,<t>,L)  =*  q^  (p,(t>,0)  =  -k^(no)outG(p,$,L)/4iri. 

We  have  introduced  the  modified  notation  (n  ) .  and  (n  )  to  call  attention 

o  m  o  out 

to  the  different  expressions  that  apply  for  in-plane  and  out-of-plane 
scattering.  Additionally,  note  that  in  the  usual  case  which  the  dimensions  of 
the  scattering  area  are  much  larger  than  those  of  the  scatterer,  both  G(4>,L) 
and  G(p,<}>,L)  can  be  replaced  by  a  single  Fourier-transform 

G(K(*))  =  /  /  dxdy  g(x,y)e1-(<^)'-  ,  (10) 


where  K ( >  =  3q(1,0,0)  -  Eo(cos<J>,sin4>,0)  for  in-plane  scattering  or 
K(<j>)  *  6o(l,0,0)  -  6(cos<fi ,sin<{> ,0)  for  out-of-plane  scattering.  That  is, 

K(<(>)  is  the  scattering  vector. 

Let  us  consider  a  scattering  area  that  consists  of  N  scattering  cen¬ 
ters  located  at  _r  =  _r..,  i  =  1,2,..N.  With  =  r- r_  ,  we  may  write 

N 

g(x,y)  =  l  Y, (q. .  q?.)  (11) 

i=1  l  ii. 

where  describes  the  scattering  perturbation  in  the  vicinity  of  the 
i  th  scattering  center.  The  widths  Aq^  and  ^2^  over  which  y^  is  significantly 
different  from  zero  describes  the  size  of  the  scattering  centers.  Inserting 
Eq.  (11)  into  Eq.  (10),  we  find 

G(K)  =  l  e1-'-!  f  S  do  y  (q.)e-^i  (12) 

^  —  OO  —00  ^ 

For  identical  scatterers  having  a  width  a_  in  the  x  direction  and  b  in  the 
y  direction,  we  find 
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-tv 

G(K)  =  l  e  - -i  ab  sine  (K  a/2ir) sine (K  b/2ir) 
i=l  y 


(13) 


CALCULATIONS  OF  SCATTERED  POWER 


The  amount  of  power  scattered  in  the  plane  of  the  waveguide  into  a 
range  of  angles  d<J>  centered  at  is  calculated  from  the  formula 

dpin(*)  =  P(2ir/0o)  c Cd>,L)c*(<|>,L)d4>  ,  (14) 

where  P  is  the  initial  power  per  unit  width  of  the  beam  in  the  waveguide. 
Similarly,  the  light  scattered  out  of  plane  into  a  range  of  angles  d<j>  centered 
at  <(>  and  into  a  range  of  radiation  modes  dp  centered  at  p  is 

d2Pout(p,<}>)  =  Prad(2iT/B)q(p,<|>,L)q  (p  ,$  ,L)dpdt|> .  (15) 

By  way  of  Eq.  (9)  part  of  the  angular  dependence  of  each  of  these  expressions 

2 

depends  on  the  factor |g|  •  For  identical,  randomly  positioned  scatterers, 

N 

|G  [  2  =  a2b2sinc2(K  a/2Tr)sinc2  (K  b/2n)  ][  e^— i  — j\  (16) 

X  y  ^=1 

For  scattering  angles  such  that  | K(<J>)  |  1  is  comparable  to  sample  dimensions, 

2 

the  sum  in  Eq.  (16)  has  the  value  N  .  Otherwise,  it  fluctuates  about  the 
average  value  N  and  we  have 

| G |  2  =  Na2b2sinc2(K  a/2w) sinc2(K  b/2ir)  .  (17) 

x  y 

From  another  point  of  view,  that  in  which  the  scattering  distribution  is  a 
random  function  g(x,y),  having  correlation  lengths  a  and  b  in  the  x  and  y 
directions,  respectively,  the  evaluation  of  | G |  is  performed  using  methods 
discussed  in  Ref.  (3).  The  result  is  similar  to  that  of  Eq.  (17),  and  is 

| G | 2  -  4  Na2b2  L(K  a)L(k  b)  (18) 

x  y 
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where 


L(q)  =  (1  +  qV1  (19) 

is  the  Lorentzian  function. 

Despite  the  differences  in  models  employed  in  the  calculations, 

2 

Eqs.  (17)  and  (18)  for  |g|  show  rather  similar  behavior.  We  will  write 

J G |  2  =  4Na2b2H((j>)  (20) 

to  separate  the  model-dependent  angular  variation  of  scattering  from  the 
dependence  on  scattering  center  size  and  number.  From  Eq.  (18),  which  we 
will  use  to  be  specific, 

H ( <}> ) -*H  ( <(> )  -  fl  +  Bo2a2(l-cos<t>)2r1  [1  +  Bo2b2sin2<}>] "1  (21) 

for  in-plane  scattering,  while 

H(<f>)->H(p ,#)  -  [1  +  (Bo-Bcos4.)2a2]_1tl  +  B2b2sin2<j>]“1  (22) 

for  out-of-plane  scattering. 


In-Plane  Scattering 

The  next  step  in  the  evaluation  of  scattered  power  is  the  calculation 
of  the  perturbation  terms  q  .  These  are  different  for  in-  and  out-of -plane 
scattering  and  for  volume  and  surface-roughness  scattering.  For  in-plane 
scattering  we  find  from  the  definitions  stated  below  Eqs.  (6)  and  (7) 

(n  ).  =  2n06n  (volume  scattering),  (23) 

o  In  i 

(n  ),  =  (n?  -  n2)<r  E2(0)  (0  /2u)p  P)  (surface  roughness  scattering). 

'o  in  l  o  o  o  o  o 
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Since  neither  of  these  expressions  depends  on  <f> ,  the  angular  distribution  of 
scattered  energy  for  in-plane  scattering  is  that  of  the  function  (<f> )  alone. 
Therefore  it  is  not  possible  to  differentiate  between  surface  and  volume 
scattering  on  the  basis  of  the  angular  distribution  of  scattered  light.  How¬ 
ever,  the  dependence  of  (n0)^n  on  the  field  at  the  surface  Eq(0)  for  scatter¬ 
ing  from  surface  roughness  shows  that  the  strength  of  scattering  from  this 
source  depends  on  the  modal  characteristics  of  the  waveguide.  This  is  not 
true  for  volume  scattering,  at  least  to  the  extent  that  the  scatterers  extend 
well  into  the  substrate, as  our  model  assumes. 

For  a  single-mode  waveguide  operating  just  above  cutoff,  an  approxi¬ 
mate  expression  for  the  square  of  the  electric  field  at  the  surface  is 

E2(0)  *  (4wu  P/S  D) (n2  -  n2)(n2  -  n2)'1  .  (24) 

o  o  o  1  l  1  o 

The  square  of  the  field  reduces  to  about  60%  of  this  value  for  a  strong 
single-mode  waveguide,  operating  just  below  the  onset  of  the  second  mode. 

By  combining  Eqs.  (9) , (14) , (20) , (21) , (23) ,  and  (24),  we  find  the  in-plane 
distribution  of  scattered  power  to  be 

[dPin(<l>)1surf  =  d<f>P(41r2/n1)(Na2b2a2/D4Ao)(n2-n2)Ho((j,)  (25) 

For  volume  scatterers,  the  corresponding  formula  is 

[dPin(,j>)]vol  =  d<fP(16,r2n2/nl)(Na2b26n2/Xo3)Ho(<!,)  *  (26) 

The  differential  in-plane  scattering  cross  section  is  defined  by  the 
expression 


ain^)  =  (dPin/d*)/PWL, 


(27) 


■ 


and  is  equal  to 

[Oin(^]surf  ■=  <4it2/ni)(Ea2bV/D4Xo)(Vn2)Ho(^  *  (28) 

[®in  C$)  ^vol  =  , 

where  Z  =»  N/LW  is  the  surface  density  of  scattering  centers.  For  scattering 
centers  of  comparable  size  and  density  in  LiNbO^,  we  have  the  ratio 

[ o .  ]  . / [a .  ]  ,  =  440(D46n2/X2a2) 

1  inJvol  in  surf  o  o 

-2  -2 

Taking  D  =  2  urn,  =  0.633  pm,  aQ  =  10  pm,  and  Sn  =  10  ,  we  find  that 

volume  scattering  centers  are  18000  times  as  effective  as  surface  roughness 
for  causing  in-plane  scattering. 


Out-of-Plane  Scattering 

The  out-of -plane  scattered-energy  distribution  is  calculated  from  the 

results  of  Eqs.  (9),  (15),  (20)  and  (22).  It  is,  however,  first  necessary  to 

obtain  expressions  for  the  perturbation  term  (hQ)  t  that  appears  in  Eq.  (9). 

The  definition  of  (n  )  from  the  discussion  below  Eqs.  (6)  and  (7)  is 

o  out 


(h  )  „.  =  (e/2wii  P  ,)I(P) 
o  out  o  rad 


g(x,y)  I(p)  = 


"ft 

/  dzE  (p,z)n(x,y,z)Eo(z) 

— oo 


(29) 


For  volume  scattering  centers,  n  =  2n26ng(x,y)  in  the  substrate  and  n  =  0  in 
the  superstrate.  The  last  part  of  Eq.  (29)  becomes 

0  * 

I(p)  1  =  +  2n2^n  /  E  (P,z)EQ(z)dz.  (30) 

—  OO 

Most  of  the  contribution  to  this  integral  comes  from  the  waveguide 
layer  -  D  £  z  <_  0.  In  this  region  Eq(z)  ^  cos(irz/D  +  tt/2)  for  a  well-confined 
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2  2  2  2  2  1  /2 

single  guided  mode,  while  E(p,z)  ^  cos[ (kQn*-kon2  +  p  )  (z  +  d/2)]  for 

even  radiation  modes.  For  large  values  of  p,  the  integrand  of  Eq.  (30) 

is  rapidly  oscillating  in  the  region  of  interest,  and  I(p)  ,  is  small.  This 

vol 

result  indicates  that  scattering  from  guided  modes  tends  to  populate  radiation 
modes  with  small  values  of  p.  These  are  modes  that  have  plane-wave  components 
that  propagate  in  the  substrate  at  small  angles  to  the  incident  beam,  and  decay 
exponentially  in  the  superstrate.  Using  the  orthogonality  between  guided  and 
radiation  modes,  Eq.  (30)  may  be  rewritten 

00  * 

1 (p ) vol  =  “2n2  5n  /  E  (P.z)EQ(z)dz  ,  (31) 

o 

where  the  integral  is  carried  out  over  the  superstrate.  For  the  case 

p«kQn2  and  n.,  =  n^,  both  guided  and  radiation  modes  have  decay  constants 

2  2  1/2 

close  to  kQ(n^-no)  ,  The  result  of  integrating  Eq.  (31)  is  then 

ivoi  =  -n26nE*(p,0)Eo(0)/ko(n^-nV/2  .  (32) 

For  scattering  from  surface  roughness,  the  perturbation  of  Eq.  (3)  is 
to  be  inserted  in  Eq.  (29),  with  the  result 

'surf  ‘  <”l  '  %)«/(p'0)Eo(0>-  (33> 

Both  expressions  are  seen  to  depend  on  the  value  of  the  fields  at 

2  (18) 

the  surface.  An  approximate  value  of  E  (0)  is  given  in  Eq.  (24).  Marcusev  ' 

ie 

derives  an  expression  for  E  (p,0)  which  we  approximate  as 
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(34) 


1 


E2(p,0)  =  (4o,MoPrad/TTg)p2/ko2(n2-n2). 


By  combining  Eqs .  (9) , (15) , (20) , (22) , (29) ,  and  (32-34),  we  obtain 
the  differential  scattered  power 


2 


(35) 


2  2  1/2 

D  =  n0<$n/k  (n,  -n  )  (volume  scattering) 
pert  2  o  1  o 


2  2 

D  =  (n,-n  )a  (surface-roughness  scattering), 
pert  loo 


The  radiation  mode  label  p  is  related  to  the  scattering  angle  9  of  Fig.  1 
according  to 


p  =  k  n~cos0 .  (36) 

o  2 

This  permits  us  to  recast  Eq.  (35)  in  the  more  common  format  of  a  differential 

2  2 

scattering  cross  section,  ^out  =  <1  P/dfi  /PWL.  We  find 

a  „  =  (32TT2n2)(j]a2b2/X4)(D2  /D)  (n2-n2)  (n2-n2)"2cos20H(0 ,<(,) ,  (37) 

out  z  o  pert  l  /  l  o 

H ( 0 , 4> )  =  [1  +  k2n2a2  (l-sin0cos<}i)2]  "'"[l  +  k2n2b2sin20sin  4> ] 

where  1  =  N/WL  is  the  surface  density  of  scattering  centers. 

For  scattering  centers  of  comparable  size  and  density  in  LiNbO^,  we 
have  the  ratio 


[a  ]  ,/[o  ]  ,  =  (0.1)  5n2/k2o2  (38; 

1  out Jvol' 1  out J surf  o  o 

-1  -2  -2 

Taking  k  *  10  yin  o  =  10  ym  and  5n  =  10  ,  we  find  that  surface  scatter- 

o  o 

ers  are  1000  times  as  effective  as  volume  scatterers  in  producing  forward, 
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out-of-plane  scattering  into  the  substrate.  This  is  opposite  to  the  situation 
that  we  found  to  exist  in  the  case  of  in-plane  scattering,  where  volume  scat¬ 
tering  centers  exerted  the  dominant  effect. 


Waveguide  Attenuation 

A  theoretical  expression  for  the  waveguide  attenuation  coefficient 
is  found  by  integrating  the  differential  scattering  cross  section  aout(0,4>) 
over  all  solid  angles  in  the  substrate.  A  graphical  analysis  of  the  integral 
suggests  that  the  dominant  contribution  comes  from  the  angular  range 

_i  /o 

-cos0  <  (k  n,a)  (39) 

'w  O  c. 

|  sin<b  I  <  (kQn2b)  1 

We  will  crudely  evaluate  the  integral  by  setting  H(0,<j>)  =  1  in  this  angular 
range  and  H(0,<(>)  =  0  elsewhere.  This  should  at  least  give  us  the  correct 
dependences  of  attenuation  on  the  parameters  kQn2a  and  k  n^b.  We  find 

“surf  =  (^/3)(2/ir)1/2n2"1/2(n2-n22)0:a1/2bo2/X^/2D)  ,  (40) 

a  =  (l/37r2)(2/7r)1/2n2/2(nJ-n2)(Za1/2bXy2/D)(n^-n2)"3 

VOl  i.  1  z  o  1  o 

for  the  attenuation  coefficients  associated  with  surface  and  volume  scattering 
centers,  respectively. 

Using  the  Theory  to  Describe  Experimental  Data 

Although  Eq.  (40)  may  be  applied  directly  to  the  analysis  of  experi¬ 
mental  attenuation  data,  Eq.  (28)  for  the  in-plane  scattered-energy  distribu¬ 
tion  must  be  extended  to  account  for  the  method  usually  employed  to  take  in¬ 
plane  scattering  data.  This  method  is  illustrated  by  Fig.  3  on  page  21. 

Light  scattered  over  a  range  of  angles  is  projected  by  a  lens  onto  a  slit  of 
width  s.  If  the  focal  length  of  the  lens  is  f  and  the  index  of  the  waveguide 
mode  is  approximately  n^,  the  range  of  angles  observed  is 
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(41) 


A<J>  =  s/n^f  , 

where  A$  is  measured  in  the  waveguide  medium.  The  scattered  power  incident 
on  the  detector  is 

AP  =  PWLo^UWn^  ( 

This  power  is  referenced  to  the  power  incident  on  the  detector  when  the  un¬ 
scattered  beam  is  projected  onto  the  slit.  This  power  is 


P  =  PW2s/fX  (43) 

max  o 

The  in-plane  scattered-energy-distribution  function  is  the  ratio 

^7  =  (LXo/n1W)oin(<!>)  (44) 

max 

For  the  case  of  scattering  by  volume  imperfections,  we  find  using  the  second 
of  Eqs .  (28), 

AP/Pmair  =  16ir2(L/W)(E«n2a2b2/X2)H  U)  (45) 

Similar  considerations  apply  to  the  analysis  of  out-of-plane  scattering  data. 

The  reader  interested  in  more  details  regarding  the  coupled  mode 
analysis  that  has  been  used  to  derive  the  results  of  this  section  may  find 
them  in  Ref.  (19)by  Marcuse,  as  well  as  in  Refs.  (3)  and  (18)  cited  earlier. 
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III.  EXPERIMENTAL  APPROACHES  TO  THE  STUDY  OF 
OPTICAL  WAVEGUIDE  SCATTERING 

An  ideal  situation  is  that  in  which  sources  of  scattering  are  ob¬ 
served  directly  in  nondestructive  fashion  with  sub-micron  resolution.  The 
electron  microscope  is  the  only  instrument  we  know  of  with  this  capability, 
and  even  this  tool  is  inadequate  when  the  scatterers  constitute  a  small  pertur¬ 
bation  of  the  average  environment.  Not  surprisingly,  this  is  the  case  with 
many  high-quality  waveguide  materials.  As  a  result,  the  major  amount  of  in¬ 
formation  that  can  be  gathered  regarding  optical  waveguide  scattering  comes 
to  us  through  scattered-energy-distribution  experiments  such  as  those  de¬ 
scribed  in  this  section.  In  these  experiments,  light  propagating  in  a  wave¬ 
guide  is  scattered  and  detected  as  a  function  of  angle  or  as  a  function  of 
the  position  at  which  scattering  occurs,  and  the  measurements  are  interpreted 
according  to  reasonable  models  to  provide  indirect  information  about  the 
nature  of  scattering  sources. 

IN-PLANE  SCATTERING  EXPERIMENTS 


Figure  3  shows  the  experimental  configuration  used  to  measure  the 
in-plane  scattered  energy  distribution  in  Ti:LiNbO^  waveguides.  An  approxi¬ 
mately  Gaussian  beam  is  prism-coupled  into  and  out  of  the  waveguide  and 
fr.cused  by  a  simple  lens  onto  a  slit  coupled  to  a  fiber  bundle.  The  slit  is 
scanned  across  the  focal  plane  and  collects  light  scattered  over  an  angular 
range  of  several  degrees  about  the  initial  beam  direction.  When  the  slit  is 
translated  a  distance  x  from  the  focus,  light  passed  through  it  to  the  detec¬ 
tor  corresponds  to  light  scattered  at  an  angle 

=  x/fng  (46) 

in  the  plane  of  the  waveguide,  where  f  is  the  focal  length  of  the  lens,  n^ 
is  the  effective  index  of  the  waveguide  mode,  and  <p  corresponds  to  the  angle 
in  the  waveguide  material. 

For  the  better  waveguides  we  have  worked  with,  the  detected  signal 

at  n  $  =  1°  is  down  40  dB  or  more  from  the  peak  signal  at  the  focus.  A 

g 
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Fig.  3.  Experimental  configuration  for  measuring  the 
in-plane  scattered-energy  distribution. 


photomultiplier  tube  is  thus  a  convenient  detector,  since  it  possesses  the 
requisite  dynamic  range.  Our  measurements  are  facilitated  by  use  of  a  lock-in 
amplifier  operating  in  a  frequency  selective  mode.  The  output  from  the  lock-in 
is  converted  to  logarithmic  format  by  a  log-amplifier  and  displayed  using  an 
x-y  recorder. 

In-plane  scattering  data  are  obtained  in  two  steps.  First,  the 
lock-in  is  employed  on  a  high -sensitivity  setting  to  record  scattered  light 
levels  in  a  range  -30  dB  to-50  dB  below  the  peak  intensity  of  the  unscattered 
focused  beam.  Then  the  sensitivity  of  the  lock-in  is  reduced  and  a  scan  is 
made  of  the  intensity  profile  of  the  unscattered  focused  beam. 

For  the  in-plane  scattering  data  presented  in  this  report,  the 
scattered  intensity  measured  in  dB  reduction  from  the  peak  intensity  is 
plotted  as  a  function  of  the  scattering  angle  measured  external  to  the  wave¬ 
guide,  <(>  =  n  <)>.  This  is  a  convenience,  since  n  varies  between  2.2  and 

S  S 

2.3  in  our  experiments,  depending  on  mode  polarization  and  direction  of  beam 
propagation.  Interpretation  of  the  data  must  also  be  guided  by  the  understand¬ 
ing  that  the  scattered  intensity  increases  linearly  with  the  ratio  L/W,  where 
L  is  the  length  of  the  scattering  medium,  and  W  is  the  beam  width.  In  our  ex¬ 
periments,  L  =  18  mm  and  W  =  1.5  ram 

In-plane  scattering  results  to  be  described  in  Sec.  IV  were  obtained 
using  the  system  of  Fig.  3,  with  one  modification;  namely,  a  mechanical  chopper 
was  used  in  place  of  the  oscillating  mirror  to  modulate  the  beam.  A  potential 
advantage  of  the  oscillating  mirror  is  that  light  scattered  at  the  input  prism 
is  not  modulated  to  the  extent  that  light  scattered  in  the  waveguide  is  and 
can  therefore  be  rejected  at  the  detector.  However,  there  is  still  the  problem 
of  modulated  scattering  at  the  output  prism.  This  constitutes  a  source  of 
noise  that  must  be  evaluated.  Early  in  the  program  we  verified  that  waveguide 
scattering  dominated  over  prism  scattering  in  the  case  of  a  fairly  poor  wave¬ 
guide.  Recent  evidence  suggests  that  waveguide  quality  is  improved  to  the 
level  where  prism  effects  can  influence  the  data.  At  this  level  it  is  appro¬ 
priate  to  take  data  with  zero  input-output  prism  separation  and  to  subtract 
the  observed  scattering  distribution  from  that  obtained  with  a  large  prism 
separation.  Also  data  should  be  taken  for  several  beam  locations  within  the 
waveguide,  and/or  several  placements  of  the  coupling  prisms,  and  the  most 
representative  data  used  for  more  detailed  analysis.  The  data  presented  in 


Sec.  IV  does  not  for  the  most  part  reflect  these  recommendations.  Justifi¬ 
cation  is  to  be  found  in  the  nature  of  the  data  and  the  methods  employed  to 
interpret  it.  First,  we  avoid  quantitative  analysis  of  results  which  could 
be  in  error  because  of  prism-associated  scattering.  Second,  we  emphasize 
changes  in  the  scattered-energy  distribution  that  result  from  changes  in  ex¬ 
perimental  conditions.  Prism  coupling  in  and  out  of  the  waveguide  is  a  con¬ 
stant  of  each  experiment,  and  therefore  does  not,  presumably,  influence  the 
changes  in  scattering  that  are  relevant  to  our  analysis. 

OUT-OF-PLANE  SCATTERING  EXPERIMENTS 


Early  in  the  program  considerable  effort  was  directed  toward  out-of¬ 
plane  scattering  measurements,  particularly  the  related  measurement  of  wave¬ 
guide  attenuation.  The  technique  that  was  employed  most  often  is  shown  in 
Fig  4.  It  was  thought  that  the  technique  offered  several  advantages  over  more 

conventional  means  for  measuring  waveguide  attenuation,  excluding  pyroelectric 
(7  8) 

means.  ’  In  any  case,  a  nonpyroelectric  method  is  needed  for  the  study  of 
the  out-of-plane  scattered- energy  distribution  in  LiNbO^. 

In  the  procedure  of  Fig.  4  an  oscillating  mirror  is  employed  to  modu¬ 
late  coupling  into  the  waveguide,  with  the  idea  that  light  scattered  in  the 
prism  will  not  be  similarly  modulated  and  can  therefore  be  rejected  at  the 
detector.  Second,  light  scattered  into  the  substrate  is  studied  with  the  idea 
that  the  scattered  intensity  in  the  substrate  will  be  greater  than  in  the 
superstrate,  thus  affording  a  signal-to-noise  advantage  important  to  the  study 
of  very-low-loss  waveguides. 

We  found  that  light  scattered  at  the  prism  was  not  totally  modula¬ 
tion  free,  and  that  much  of  it  was  directed  into  the  substrate  where  it  was 
necessary  to  shield  it  from  the  fiber  bundle.  Shielding  was  also  necessary 
at  the  output  end  of  the  waveguide,  whether  or  not  an  output  coupling  prism 
was  employed.  The  amount  of  shielding  required  increases  with  the  substrate 
thicknesses  and  refractive  index.  Thickness  is  a  factor  because  it  consti¬ 
tutes  the  minimum  distance  from  the  prism  that  shielding  can  be  positioned. 
Refractive  index  is  a  factor  because  it  influences  the  number  of  reflections 
that  light  undergoes  before  leaving  the  substrate  and  being  blocked  by  the 
shielding. 
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Fig.  4.  Experimental  configuration  for  observing 
out-of-plane  waveguide  scattering. 


The  need  for  shielding  reduces  the  overall  length  of  the  waveguide 
that  can  be  scanned  to  measure  attenuation.  This  reduces  any  advantage  that 
may  be  associated  with  increased  scattered  intensity  in  the  substrate  region. 
Consequently,  a  geometry  in  which  light  scattered  into  the  superstrate  is 
detected  may,  in  fact,  be  superior  to  the  geometry  shown  in  Fig.  4.  For  the 
measurement  of  waveguide  attenuation.  The  geometry  of  Fig.  4,  however,  would 
still  be  relevant  for  the  study  of  the  out-of-plane  scattered  energy  distribu¬ 
tion. 

We  suspended  our  measurements  in  out-of-plane  scattering  in  favor  of 
in-plane  measurements  that  were  straightforward  to  perform  and  still  provided 
us  with  a  wealth  of  information.  Later  in  the  program,  theoretical  calcula¬ 
tions  showed  that  sources  for  out-of-plane  and  in-plane  scattering  could  be 
very  different.  In  view  of  the  relevance  of  in-plane  scattering  to  the 
spectrum-analyzer  application,  we  continued  to  focus  our  attention  on  sources 
of  in-plane  scattering  for  the  duration  of  the  program. 


IV.  IN-PLANE  SCATTERING  IN  Ti:LiNb03  WAVEGUIDES 

In  this  section  we  present  and  interpret  in-plane  scattering  data 
for  various  TirLiNbO^  waveguides.  The  data  were  obtained  using  the  experi¬ 
mental  method  described  in  Sec.  III.  TitLiNbO^  waveguides  were  emphasized  for 

study  because  they  are  the  current  waveguides  of  choice  for  the  development  of 

(1  9) 

an  integrated  optical  spectrum  analyzer.  *  The  data  to  be  presented  have 
led  us  to  the  following  conclusions  regarding  the  nature  of  scattering  in 
these  waveguides: 

(1)  Bulk  scattering,  at  least  some  of  which  is  diffusion- induced, 
is  far  more  important  than  surface  scattering. 

(2)  The  bulk  scatterers  are  anisotropic  in  refractive  index  and 
possibly  in  shape. 

(3)  The  scatterers  are  associated  with  regions  that  have  different 
electrooptic  properties  from  those  of  the  crystal  as  a  whole. 

(4)  Scattering  increases  in  proportion  to  the  amount  of  Ti  used  to 
make  the  waveguide,  but  is  not,  within  limits,  sensitive  to  the  time  and  dura¬ 
tion  of  the  diffusion  treatment. 

(5)  Scattering  is  similar  for  x-cut  and  y-cut  crystals. 

The  recipe  that  we  employ  for  our  best  waveguides  requires  the  use 

O  o 

of  thin  150  A-300  A  Ti  layers  diffused  from  3-6  h  at  950°C  in  flowing  0^,  fol¬ 
lowed  by  a  rapid  quench  to  600°C  and  a  slow  cool  to  room  temperature.  For 
waveguides  formed  using  thicker  Ti  films,  some  improvement  in  quality  is 
always  obtained  by  lightly  polishing  the  finished  waveguide,  even  when  all  Ti 
has  been  diffused  into  the  surface. 

The  above  recipe  applies  for  waveguides  used  with  beam  propagation 
perpendicular  to  the  c-axis.  For  the  parallel  configuration,  however,  our 
best  low-scattering  performance  was  obtained  using  a  thick  Ti  film  (about 

O 

720  A)  and  a  long  post-fabrication  polish — long  enough,  in  fact,  to  remove 
most  of  the  initially  formed  waveguide  layer. 


26 


IN-PLANE  SCATTERING  IN  A  LiNbCh  WAVEGUIDE 
FORMED  USING  A  280  A  Ti  FILM 


The  preferred  configuration  for  an  RF  spectrum  analyzer  employs  TE 
polarized  light  propagating  nearly  perpendicular  to  the  LiNbO^  optic  axis.  For 
this  configuration  we  have  found  the  least  scattering  in  waveguides  formed 
using  relatively  thin  Ti  films.  The  data  presented  here  are  for  a  sample  ob- 

O 

tained  by  diffusing  a  280  A-thick  Ti  layer  into  the  substrate  for  3  h  at 
950°C  in  flowing  0~.  However,  we  have  tentative  evidence  that  some  small  im- 

^  O 

provement  may  be  obtained  by  going  to  even  thinner  (150  A)  Ti  films.  The 
sample  we  will  describe  was  lightly  polished  following  fabrication  to  improve 
its  surface  condition.  This  procedure  and  its  benefits  will  be  discussed  in 
detail  later. 

Figure  5  shows  typical  in-plane  scattered-energy  distributions 
for  the  TEq  mode  for  four  different  directions  of  beam  propagation:  0° ,  90°, 
+45°,  and  -45°  degrees  to  the  optic  axis.  The  anisotropy  of  the  observed 
scattering  strongly  suggests  that  the  scattering  mechanism  is  not  surface 
roughness.  (The  surface  polish  of  the  initial  substrate  was  an  "optical" 
polish  provided  by  the  supplier  using  proprietary  methods) .  The  anisotropy 
could  result  from  volume  inhomogeneities  that,  as  seen  by  the  waveguided  beam, 
are  not  simply  circular  in  shape.  Or  the  anisotropy  could  result  from  circu¬ 
lar  scattering  centers  having  an  anisotropic  refractive  index.  Alternatively, 
both  geometric  and  refractive-index  anisotropies  could  be  relevant. 

Figure  6  shows  the  scattered  energy  distribution  for  the  same 

waveguide  when  the  TMq  mode  is  propagated  at  the  same  four  angles  to  the 

optical  axis.  The  variation  with  angle  is  similar  to  that  for  TEq  modes  in 

the  case  of  propagation  parallel  and  perpendicular  to  the  optic  axis.  However, 

TM  mode  scattering  is  least  at  +  45°  to  the  optic  axis,  while  TE  mode  scatter- 
o  —  o 

ing  is  greatest  in  these  directions. 

This  behavior  is  difficult  to  explain  if  the  scattering  centers  are 
assumed  to  have  an  isotropic  refractive  index  and  an  anisotropic  shape.  For 
example,  rod-like  scatterers  with  their  long  axis  at  90°  to  c_  could  be  pro¬ 
posed  to  explain  the  data  taken  for  propagation  parallel  and  perpendicular  to 
the  optic  axis.  However,  scattering  at  +  45°  would  then  be  expected  to  be 
intermediate  to  that  at  0°  and  90°,  in  contrast  to  what  we  observe. 
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Fig.  5.  In-plane  scattered  energy  distributions  for  four 
different  directions  of  TEq  mode  propagation. 

(NOTE:  The  position  of  the  arrowhead  of  the  c-axis  was  not 
deemed  to  be  experimentally  significant,  and  could 
be  opposite  that  shown  by  the  insets) . 
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Fig.  6.  In-plane  scattered  energy  distributions  for  four 
different  directions  of  TMq  mode  propagation. 

(NOTE:  The  position  of  the  arrowhead  of  the  c-axis  was  not 
deemed  to  be  experimentally  significant,  and  could 
be  opposite  that  shown  by  the  insets) . 
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More  likely  is  the  possibility  that  the  scattering  centers  exhibit 
an  anisotropic  refractive  index  variation.  Regions  where  the  crystal  stoichiom¬ 
etry  fluctuates  or  where  crystalline  compounds  containing  Ti  are  formed  during 
diffusion  are  possible  examples  of  such  anisotropic  scattering  centers.  The 
simple  process  of  Ti  diffusion  is  known  to  produce  anisotropic  index  changes, 
even  in  the  absence  of  compound  formation,  because  of  the  tendency  of  the  Ti 
to  occupy  well  defined  lattice  sites. 

IN-PLANE  SCATTERING  IN  A  LiNbOi  WAVEGUIDE 
FORMED  USING  A  720  X  Ti  FILM 


The  Surface  Contamination  Layer 

We  have  observed  considerably  greater  scattering  than  seen  in  Figs. 

5  and  6  in  the  case  of  LiNbO^  waveguides  formed  using  thicker  Ti  films.  A 
quantitative  indication  of  the  variation  of  waveguide  quality  with  Ti-film 
thickness  is  shown  by  the  graph  of  Fig.  7.  Increased  scattering  is  associ¬ 
ated  with  the  increase  in  a  granular  "orange  peel"  layer  found  on  the  surface 
of  the  sample  after  Ti  diffusion.  Figure  8  shows  an  example  of  this  granular 
layer,  taken  using  differential  interference  contrast  (Nomarski)  microscopy. 
The  demarkation  between  the  granular  diffused  region  and  an  adjacent  undiffused 
region  is  readily  apparent.  Observation  of  the  demarkation  using  an  inter¬ 
ference  microscope  suggests  that  the  thickness  of  the  granular  layer  is  a  few 
hundred  Angstroms;  i.e.,  comparable  to  the  thickness  of  the  diffused  metal  film. 
The  characteristic  lateral  dimension  of  the  granular  structure  is  about  one 
micron.  The  features  seen  in  Fig.  8  are  roughly  characteristic  of  all  wave- 

O 

guides  we  have  formed  using  a  Ti-film  thickness  of  720  A.  For  thinner  films, 
the  granular  structure  is  much  less  apparent,  but  still  present.  Fig.  9 
shows  the  surface  texture  obtained  in  an  extreme  case,  that  of  a  26  h  diffu- 

a 

sion  of  a  1000  A-thick  Ti  film.  The  granular  structure  is  approximately  two 
microns  in  size.  We  emphasize  that  the  structure  is  not  undiffused  Ti02,  but 
rather  is  associated  with  the  diffused  species. 
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Nomarski  micrograph  showing  the  granular 
surface  texture  associated  with  Ti  dif¬ 
fusion  (750  A  Ti  film,  3  h  diffusion, 

500X  magnification).  Vertical  line  shows 
the  demarkation  between  the  diffused  region 
(right)  and  the  undiffused  region  (left) . 


Fig.  9.  Nomarski  micrograph  showing  the  granular 
surface  texture  associated  with  Ti  dif¬ 
fusion  (1000  A  Ti  film,  26  h  diffusion, 

500X  magnification).  Vertical  line  shows 
the  demarkation  between  the  diffused  region 
(right)  and  the  undiffused  region  (left). 
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Waveguide  Improvement  by  Post-Diffusion  Polishing 

It  is  not  surprising  that  diffusion  of  the  metal  increases  the 
volume  of  the  crystal;  however,  we  can  only  speculate  as  to  why  a  nonuniform 
surface  profile  results.  Burns  et  al  ^^have  obtained  evidence  of  higher 
than  anticipated  Ti  concentration  within  the  top  0.3  im  of  diffused  LiNbO^ 
waveguides  having  an  overall  depth  of  several  microns.  They  tentatively 
ascribe  this  to  the  formation  of  Li-Ti-0  compounds  such  as  L^Ti^O^.  Per¬ 
haps  the  granularity  we  observe  is  evidence  of  the  formation  of  crystallites 
of  this  material.  The  crystallites  would  presumably  have  an  index  consider¬ 
ably  different  from  that  of  LiNbO^,  and  so  would  be  efficient  volume  scatter¬ 
ing  centers.  In  corroboration  of  this  expectation  we  find  sometimes  consider¬ 
able  reduction  in  waveguide  scattering  by  polishing  the  waveguide  surface. 

The  polishing  is  performed  after  diffusion  using  a  colloidal  suspension  of 
silica  (trade  name:  Syton)  in  conjunction  with  a  microcloth  lap.  The  polish¬ 
ing  is  done  by  applying  several  pounds  per  square  inch  pressure  to  the  sample 
as  it  moves  about  the  circumference  of  a  vibrating  table.  Polishing  continues 
for  at  least  a  time  long  enough  to  remove  the  granular  "orange  peel"  structure 
on  the  surface.  In  our  experiments,  this  was  45-90  minutes.  Polishing  for 
this  length  of  time  does  not  appear  to  significantly  alter  the  modal  properties 
of  the  waveguide. 

We  have  found  continued  scattering  reduction  in  waveguides  formed 

O 

using  720  A  Ti  films  when  we  continue  to  polish  the  waveguide  surface  for 
times  well  in  excess  of  those  required  simply  to  remove  the  granularity.  How- 

O 

ever,  for  waveguides  formed  using  280  A  Ti  films,  there  appears  to  be  no  advant- 

O 

age  to  extra  polishing  once  the  granularity  is  removed.  And  for  150  A  films, 
where  the  granularity  is  barely  visible,  there  may  be  no  advantage  to  any 
post-diffusion  polishing. 

Figure  10  and  11  show  the  dramatic  improvement  in  waveguide 

quality  that  can  result  from  long  polishing  of  samples  formed  using  thick  Ti 

films.  The  waveguide  employed  in  the  experiment  was  fabricated  by  diffusing 
0 

720  A  of  Ti  into  the  y-cut  surface  for  3  h  at  950°C  in  flowing  O2,  followed 
by  a  rapid  quench  to  600°C.  The  data  show  the  typical  peak  scattering  level 
in  the  waveguide,  measured  at  an  external  angle  of  1°,  as  a  function  of  total 
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polishing  time.  Data  for  the  case  of  propagation  perpendicular  and  parallel 
to  the  optic  axis  are  presented  in  Figs.  10  and  11  respectively.  Through 
7  h  polishing,  there  is  monotonic  reduction  in  scattering  of  TM  modes,  while 
TE  modes  remain  essentially  static.  Perhaps  this  reflects  a  greater  sensi¬ 
tivity  of  TM  modes  to  surface  conditions,  owing  to  the  larger  value  of  the 
electric  field  at  the  surface  for  these  modes. 

After  7  h  polishing,  the  scattering  levels  for  the  three  modes  that 
propagate  as  ordinary  waves  (TEq,  k  | )  c;  TMq,  k  | |  £;  and  TMq,  k  J_  c)  are 
comparable  and  significantly  less  than  the  scattering  level  for  the  one  mode 
that  propagates  as  an  extraordinary  wave  (TEq,  k  J_  c) .  This  suggests  that 
the  scatterers  could  have  an  anisotropic  refractive  index,  and  share  the 
same  symmetry  axis  as  the  host  crystal.  It  is  also  requisite  that  the  scat¬ 
ters  be  roughly  isotropic  in  shape;  that  is,  roughly  circular  when  viewed  in 
the  plane  of  the  waveguide,  if  this  model  is  to  apply. 

After  11  h  total  polishing,  a  dramatic  reduction  in  scattering  level 
of  TEo  modes  is  shown  by  the  data.  (Unfortunately,  TMq  modes  were  not  examined). 
Scattering  for  the  case  of  propagation  parallel  to  the  optic  axis  remains  10  dB 
lower  than  scattering  for  the  case  of  propagation  perpendicular  to  the  optic 
axis,  but  both  experiments  show  6-7  dB  reduction  in  scattering  level  resulting 
from  the  additional  4  h  polishing.  A  possible  interpretation  is  that  the  Ti- 
rich  surface  contamination  layer  suggested  by  the  data  of  Burns  et  al  has 
been  completely  polished  away,  leaving  a  superior  waveguide. 

Figure  12  shows  the  scattered-energy  distribution  for  the  case 
TEq,  k  | |  c_,  11  h  polish,  which  is  the  best  waveguide  scattering  performance 
we  have  obtained  to  date.  The  performance  of  the  same  waveguide  when  used 
for  propagation  perpendicular  to  the  optic  axis,  however,  was  no  better  than 

O 

that  shown  in  Fig.  5,  for  a  waveguide  which  was  formed  with  a  280  A  Ti  film, 
and  which  did  not  have  extensive  post-diffusion  polishing. 
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Fig.  12.  In-plane  scattered-energy  distribution 
for  the  TE0  mode  propagating  parallel 
to  the  optic  axis,  after  11  h  total 
polishing  time. 
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V.  POSSIBLE  SOURCES  OF  RESIDUAL  SCATTERING  IN 
Ti:LiNb03  WAVEGUIDES 

As  Sec.  iv  has  shown,  near  surface  scattering  tentatively  identified 
with  Li-Ti-0  compound  formation  is  a  potentially  significant  source  of  scatter¬ 
ing  in  TirLiNbO^  waveguides .  Nevertheless,  it  is  not  the  only  source  of 

scattering.  We  have  considered  several  possible  sources  for  the  residual  scat¬ 
tering  in  TirLiNbO^  waveguides  that  is  not  eliminated  through  the  use  of  low 

Ti  concentrations  or  post-diffusion  polishing.  These  include  (1)  localized 

(12) 

occurrence  of  the  phase  transformation  LiNbOy+LiNb^Og;  (2)  Ti-free  regions 

resulting  from  local  imperfections  in  the  diffused  film;  and  (3)  occurrence  of 

microscopic  antiphase  domains  (microdomains) ,  which  are  regions  of  crystalline 

(13) 

imperfection  characterized  by  c^-axis  reversal.  A  fourth  possible  source 

of  the  observed  anisotropic  scattering  in  LiNbO^  is  anisotropic  surface  rough¬ 
ness  associated  with  the  uniaxial  crystalline  nature  of  the  material.  While 
this  is  a  candidate  worthy  of  study  in  the  future,  we  have,  during  the  program, 
considered  only  the  first  three  candidate  scattering  mechanisms,  and  these  to 
varying  degrees.  We  present  our  results  in  this  section.  Emphasis  is  placed 
on  the  possible  occurrence  of  antiphase  microdomains,  since  this  scattering 
mechanism  was  given  the  greatest  attention  and  provided  the  most  intriguing 
experimental  results. 

ANTIPHASE  MICRODOMAINS  AND  ELECTRIC-FIELD 
INFLUENCE  ON  SCATTERING 

(13) 

Etching  studies  by  Ohnishi  have  shown  that  domain  reversal  can 
occur  at  the  surface  of  LiNbO^  upon  heating  at  the  high  temperatures  required 
to  make  waveguides.  In  unheated  samples,  reversed  domains  were  found  to  be 
needle-like  structures  having  a  diameter  of  less  than  1  ym  and  a  length  typi¬ 
cally  300  um,  lying  parallel  to  the  c-axis.  We  have  found  similar  effects 
upon  etching  LiNbO^  samples,  however,  the  observation  of  propagation-direction- 
dependent  scattering  in  LiNbO^  waveguides  first  led  us  to  suspect  antiphase 
microdomains  as  a  scattering  mechanism,  owing  to  their  cylindrical  shape. 

Later  we  considered  that  the  cylindrical  structures  seen  after  etching  were 
oriented  in  the  wrong  direction  to  result  in  increased  scattering  for 


d 
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propagation  perpendicular  to  the  c-axis,  as  we  had  observed.  Nevertheless,  by 
the  time  of  this  consideration  we  had  obtained  positive  indication  of  an 
electric-field  dependent  component  of  scattering  that  increased  the  attractive¬ 
ness  of  the  microdomain  hypothesis. 

Electric-Field  Dependent  Scattering  in 

Ti:LiNbO„  Waveguide 
- j - 

The  experiment  that  we  performed  is  indicated  schematically  in 
Fig.  13a.  Silver-paste  electrodes  were  applied  to  the  surface  of  a  good 
quality  waveguide,  about  1-mm  apart.  Light  was  prism-coupled  into  and  out  of 
the  waveguide  and  propagated  in  the  gap  region.  The  in-plane  scattering  dis¬ 
tribution  was  measured  with  500  V  and  0  V  applied.  The  results  of  the  experi¬ 
ment  are  shown  in  Figs.  14  and  15.  An  electric  field  applied  parallel  to  the 
optic  axis  is  seen  to  result  in  a  significant  increase  in  waveguide  scattering 
(20  dB  at  1°).  The  scattering  decreases  when  the  applied  field  is  removed. 

The  magnitude  of  the  effect  appears  to  be  independent  of  the  polarity  of  the 
applied  voltage. 

Figure  13b  shows  how  this  observation  may  be  interpreted  using  the 

microdomain  hypothesis.  In  the  absence  of  an  applied  field,  the  index  of  the 

microdomain  is  close  to  that  of  its  surroundings  and  minimum  scattering  occurs. 

In  the  presence  of  an  applied  field  an  electrooptic  index  change 
3 

An  =  0.5  ner33EQ  is  induced  in  the  crystal.  Owing  to  c-axis  reversal  in  the 

microdomain,  the  sign  of  the  index  change  in  that  region  is  opposite  to  that 

in  the  host  crystal.  The  optical  beam  therefore  sees  a  net  index  change  of 

| 6n |  =  2 | An |  localized  in  the  vicinity  of  the  microdomain.  The  sign  of  6n 

depends  on  the  polarity  of  the  applied  voltage;  however,  scattering  from  the 

2 

induced  index  change  goes  as  Sn  and  is  polarity  independent. 

Photorefractive  Effects  on  Waveguide  Scattering 

As  a  further  test  of  the  microdomain  hypothesis,  we  attempted  to 
generate  a  similar  scattering  effect  to  that  shown  in  Fig.  14  using  the  photo- 
refractive  properties  of  LiNbO^.  In  this  experiment,  a  beam  of  wavelength 
515  nm  was  propagated  for  a  period  of  time  in  TiiLiNbO^  sample  109.  The 
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Fig.  13(b) .  Microdomain  model  for  voltage-induced 
scattering. 


exposure  was  sufficient  to  cause  significant  photorefraction.  By  this  process 

negative  space  charge  migrates  to  the  +c  side  of  the  optical  beam,  leaving 

(14) 

positive  space  charge  at  the  -c  side.  The  electric  field  associated  with 

this  space  charge  produces  a  refractive  index  change  by  means  of  the  electro¬ 
optic  effect. 

In  the  experiment  illustrated  in  Fig.  16,  exposure  was  carried  out 
with  the  515-nm  beam  propagating  parallel  to  the  optic  axis.  This  direction 
was  chosen  to  inhibit  the  formation  of  hologram  gratings  caused  by  the  inter¬ 
ference  of  scattered  and  unscattered  components  of  the  exposing  beam.  These 

(15' 

holograms  cause  additional  scattering  that  obscures  the  effects  of  interest.  " 
,By  propagating  the  exposing  beam  parallel  to  the  optic  axis,  no  additional 
scattering  was  observed. 

It  was  anticipated,  however,  that  during  exposure  space  charge  would 

develop  on  either  pole  of  the  microdomain  scattering  centers,  as  illustrated 

in  Fig.  16.  When  the  waveguide  is  subsequently  interrogated  using  a  beam 

propagating  perpendicular  to  the  optic  axis,  the  index  change  produced  by  the 

space  charge  should  result  in  additional  scattering.  The  extent  of  additional 

scattering  is  expected  to  be  on  the  order  of  that  produced  by  applied  fields 

approaching  100  V/mm,  since  it  is  known  that  space-charge  fields  of  that  magni- 

(14) 

tude  are  common  using  photorefraction.  We  therefore  expected  to  see  an 

additional  scattering  comparable  to  that  observed  in  Fig.  14.  However,  no 
discernible  change  in  scattering  resulted. 

It  appears  that  the  scattering  centers  are  electrooptically  different 
but  photorefractively  similar  to  the  host  crystal,  if  this  null  result  is  to 
be  explained.  This  is  damaging  to  the  microdomain  hypothesis.  The  hypothesis 
could  be  salvaged,  however,  if  it  were  found  that  the  space  charge  decayed  dur¬ 
ing  the  10  to  20  minutes  that  elapsed  between  exposure  and  interrogation.  This 
is  not  seen  to  be  likely,  in  view  of  the  fact  that  space-charge  holograms  re¬ 
corded  in  LiNbO^  samples  seldom  decay  this  rapidly.  Moreover  additional  scat¬ 
tering  was  not  observed  during  exposure  as  might  have  been  expected  if  there 
were  a  significant  space-charge  build  up,  even  a  temporary  one. 
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Fig.  16.  Photorefractive  Test  of  the  Microdomain 
Model 
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Other  Electrical  Scattering  Experiments  in  LiNbO^ 


We  perceive  the  electric-field  dependence  of  scattering  to  be  poten¬ 
tially  useful  as  a  means  for  identifying  the  source  of  residual  scattering  in 
LiNbO^,  despite  the  above  described  set  back  to  the  microdomain  hypothesis. 
Accordingly  we  performed  several  other  experiments  of  an  electrical  nature  with 
the  hope  of  understanding  the  mechanism  at  work.  The  first  experiment  was  an 
attempt  to  see  the  electrically  induced  scattering  using  Nomarski  microscopy. 
This  was  unsuccessful,  even  after  we  modified  our  Zeiss  instrument  for  use  with 

transmitted  light.  The  reason  is  to  be  found  in  the  small  refractive-index 
-4 

change  (=  10  )  and  in  the  small  path  length  (=  1  urn)  over  which  the  phase  per¬ 

turbation  must  accumulate  when  viewed  normal  to  the  waveguide  surface. 

Then  we  examined  scattering  in  a  bulk  sample  of  LiNbO^  about  3-mm 
thick  and  2.3  mm  wide.  This  enabled  us  to  avoid  waveguide  mode  effects  and 
f ringing-field  effects  associated  with  surface  electrodes.  We  measured  scat¬ 
tering  at  522  V/rnm  and  0  V/mra  and  observed  no  difference.  Later  we  subjected 
the  sample  to  an  anneal  treatment  at  950°C  with  the  idea  that  scattering  cen¬ 
ters  might  be  generated  at  the  high  temperatures  employed  for  waveguide  fab¬ 
rication.  Yet,  we  continued  to  observe  no  voltage  induced  scattering  in  this 
sample . 

Next,  we  subjected  the  sample  to  an  1.5  h  anneal  at  1140°C  to 
purposely  introduce  microdomains  by  exceeding  or  at  least  approaching  the 
Curie  temperature.  As  a  result  of  the  heat  treatment  the  sample  acquired  a 
frosted  appearance,  indicating  that  we  achieved  our  goal.  We  attempted  to  see 
electric-field  effects  in  the  sample  using  Nomarski  microscopy.  Again  we  met 
with  no  success,  despite  the  apparent  proliferation  of  microdomains.  Then  we 
observed  scattering  in  the  sample  both  with  and  without  an  applied  voltage  of 
522  V/mm.  In  either  case  scattering  was  severe.  However,  no  voltage-depen¬ 
dent  perturbation  of  background  scattering  was  observed. 

A  final  attempt  to  test  the  microdomain  hypothesis  involved  the  ex¬ 
posure  of  an  out-diffused  waveguide  to  an  Intense  beam  of  515  nm  radiation. 

The  sample  was  photorefractively  damaged  with  the  idea  of  generating  surface 
charge  localized  near  scattering  centers.  Then  the  surface  was  viewed  with 
an  SEM  in  the  hope  of  seeing  the  effects  of  surface  fields  on  beam  deflection. 
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In  this  experiment,  the  sample  was  not  gold-coated  since  this  would  have 
shorted  the  electric  fields  of  interest.  Unfortunately  but  not  surprisingly, 
the  SEM  beam  itself  generated  surface  space  charge  that  obscured  any  interest¬ 
ing  effects.  The  space  charge  accumulated  preferentially  over  certain  regions 
of  the  crystal  surface,  probably  in  response  to  spatial  variations  in  cleanli¬ 
ness.  The  photograph  of  Fig.  17  resulted.  The  background  speckle  pattern 
could  possibly  be  associated  with  waveguide  scattering  centers;  however,  the 
position  of  the  individual  speckles  changed  with  each  new  raster  scan  of  the 
electron  beam,  indicating  otherwise. 

Our  results  suggest  that  the  voltage-dependent  scattering  in 
TiiLiNbO^  waveguides  is  associated  with  Ti-dif fusion,  and  not  with  the  high 
temperatures  of  fabrication.  They  also  reinforce  the  conclusion  to  which  we 
were  drawn  by  our  first  photorefractive  experiment,  described  earlier:  micro¬ 
domains  are  not  the  source  of  the  electric-field  dependent  scattering  in 
Ti:LiNbO^  waveguides.  Unfortunately,  the  experiments  have  not  been  suffi¬ 
ciently  informative  to  suggest  other  mechanisms  for  the  observed  effect.  Nor 
have  they  been  sufficiently  complete  to  definitely  rule  out  the  possibility  of 
microdomain  contributions  to  scattering.  Other  interpretations  of  our  present 
results  could  come  to  the  forefront  in  the  future  and  make  necessary  the  resur¬ 
rection  of  microdomains  as  a  scattering  mechanism.  However,  this  remains  to 
be  seen. 

Ti-FREE  REGIONS  AND  SCATTERING  CENTERS 
ASSOCIATED  WITH  SAMPLE  CLEANLINESS 

Owing  to  imperfect  sample  cleaning  and  the  inevitable  accumulation  of 
surface  dust  prior  to  evaporation  of  the  Ti  film  used  to  make  a  waveguide,  it 
is  likely  that  certain  regions  of  the  film  will  be  destroyed  during  the  high 
temperature  waveguide-fabrication  process.  This  results  in  Ti-free  regions  of 
the  waveguide  having  an  index  lower  than  that  of  the  diffused  regions,  and 
therefore  able  to  scatter  light. 

Ignoring  the  complications  introduced  by  Li»0  outdiffusion  during 
waveguide  fabrication,  the  index  perturbation  associated  with  the  Ti-free 
region  is  -Ang  or  -AnQ,  where  Ane(AnQ)  is  the  index  change  produced  in  the 
crystal  as  a  whole  for  extraordinary  (ordinary)  polarized  light.  The 


is  consistent 


expectation  that  Ang  is  typically  50%  larger  than  An 
with  the  trend  to  lower  scattering  levels  for  ordinary  waves,  seen  in  Figs. 

10  and  11,  but  not  consistent  with  the  magnitude  of  the  observed  difference  in 
scattering  between  ordinary  and  extraordinary  waves. 

Figure  18  shows  a  Nomarski  micrograph  of  an  edge  of  the  waveguide 
used  to  obtain  the  data  of  Figs.  5  and  6.  The  artifacts  presumably  result 
from  poor  cleanliness  at  the  sample  edge.  The  center  of  the  waveguide  is 
virtually  defect  free.  The  micrograph  reveals  circular  structures  of  5-10  pm 
diameter  that  could  be  areas  where  the  Ti-film  exploded  away  at  the  high  temp¬ 
erature.  The  structures  could  then  correspond  to  Ti-free  regions.  Also  noted 
are  smaller,  reflective  structures  which  could  correspond  to  regions  where, 
for  unknown  reasons,  Ti  was  not  exploded  away,  yet  was  not  diffused  into  the 
surface. 

A  theoretical  indication  of  the  potential  of  Ti-free  waveguide 
regions  for  scattering  is  afforded  by  the  theoretical  analysis  of  Sec.  II.  The 
relevant  formula  is  obtained  from  Eqs.  (45)  and  (21): 

AP/Pmax  =  4a/W)6n2Ia2/*2  »  (47) 

where,  in  the  notation  of  Sec.  II,  AP/P  is  the  exten*-  of  scatter  reduction 

max 

measured  at  an  external  scattering  angle  assumed  large  in  comparison  to 

X0/2na;  I  is  the  number  density  of  the  Ti-free  regions,  and  fin  is  their  average 

refractive  index  change;  L  is  the  length  of  the  scattering  medium,  and  W  is 

the  optical  beam  width.  In  our  best  waveguides  we  find  AP/P  =  10_<* 

max 

(i.e.,  -40  dB)  at  tp  =  1/60,  when  L  =  18  mm  and  W  =  1.5  mm.  Thus  we  find 
-9  2  2 

10  =  fin  Za  (48) 

_2 

Taking  a  =  10  pm  and  p  =  10  mm  (for  a  reasonably  clean  surface  by  the 
standards  of  the  program)  we  find  fin  *  0.003.  This  is  comparable  to  the  index 

change  associated  with  Ti-diffusion  in  LiNbO^,  so  indeed  the  mechanism  appears 
capable  of  causing  the  observed  levels  of  residual  scattering  in  this  material. 
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LiNbOr>LiNbiOa  PHASE  TRANSFORMATION 
AS  A  SCATTERING  MECHANISM 

Q2) 

It  is  known  that  LiNbOg  becomes  unstable  in  the  vicinity  of  850°C, 
tending  to  form  LiNbgOg  according  to  the  reaction 

3LiNb03  ->  Li20  +  LiNbgOg 

The  latter  compound  is  a  monoclinic  crystal  having  the  three  refractive  in¬ 
dices  2.28,  2.36,  and  2.40.  Islands  of  this  material  formed  during  waveguide 
fabrication  could  be  very  efficient  scatterers  of  light,  owing  to  the  signifi¬ 
cant  deviation  of  the  refractive  indices  from  the  2.2  and  2.3  values  of  LiNb0o. 
These  scatterers  are  naturally  anisotropic  and  would  also  be  likely  to  display 
electric-field  dependent  scattering  effects.  However  these  effects  would  very 
likely  be  small  in  comparison  to  the  base  scattering  level  associated  with  the 
zero-field  refractive  index  changes.  This  would  be  in  contrast  to  what  we  have 
actually  observed. 

Diffused  waveguides  in  LiNbOg  are  usually  formed  at  temperatures  of 
950°C  or  higher.  While  phase  transformation  is  not  expected  to  be  a  problem 
at  these  temperatures,  heating  and  cooling  of  the  crystal  through  the  750- 
850°C  region  could  produce  sufficient  quantities  of  LiNbgOg  to  cause  observ¬ 
able  scattering.  To  test  this  hypothesis,  a  substrate  of  LiNbOg  was  heated 
in  flowing  02  at  850°C  for  10  min.  The  sample  was  brought  into  the  hot  zone 
from  a  temperature  of  about  600°C,  and  removed  to  the  same  temperature  after 
the  heat  treatment.  Heating  to  and  cooling  from  850°C  required  several 
minutes  each  way. 

The  initially  cleaned  and  polished  sample  was  observed  to  have 
numerous  pock  marks  and  scratches  following  the  experiment.  These  were  con¬ 
veniently  observed  with  a  Nomarski  microscope  under  a  power  of  500X  (Fig.  19) . 
The  scratch  marks  ran  predominantly  perpendicular  to  the  c-axis  and  parallel 
to  the  direction  in  which  the  crystal  was  wiped  during  cleaning.  They  could, 
therefore,  reflect  strains  introduced  during  cleaning  and  polishing,  along 
which  phase  transformation  preferentially  occurred.  The  pock  marks  occurred 
in  apparent  random  locations  and  seem  to  be  the  result  of  an  explosive  re¬ 
moval  of  material  from  the  surface. 
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Fig.  19.  Nomarski  micrograph  showing  the  results  of 
a  10  min  heat  treatment  of  LiNbOj  at  850*0 
in  flowing  02*  (500X) 
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Following  this  experiment  we  examined  an  outdiffused  LiNbO^  wave¬ 
guide  that  was  made  previously  and  known  to  be  of  good  quality.  This  wave¬ 
guide  was  rapidly  quenched  from  a  temperature  of  above  950°C  to  about  600- 
700°C  as  part  of  its  fabrication.  The  length  of  time  spent  in  the  range 
750-850°C  is  not  known;  however,  it  was  as  short  as  conveniently  possible  with¬ 
out  risking  fracture  of  the  sample  owing  to  an  overly  rapid  quench.  Under 
Nomarski  observation  the  waveguide  surface  showed  none  of  the  scratches  and 
pockmarks  observed  in  Fig.  19.  Consequently,  the  problems  seen  there  are  not 
characteristic  of  good  waveguides  of  LiNbO^.  However,  an  incipient  stage  of 
these  problems  is  likely  to  be  present  and  may  possibly  be  causing  observed 
scattering  levels  in  LiNbO^- 

Note  added  during  preparation  of  final  draft;  In  recent  experiments 
we  have  reexamined  the  electric-field  dependence  of  scattering  using  electrode 
gaps  of  2  and  3  mm,  as  well  as  1  mm.  We  find  that  the  observed  electric-field 
effects  vanish  as  the  gap  becomes  wider  than  the  waveguided  beam.  This  leads 
us  to  believe  that  the  observed  effects  result  from  electric-field  nonuniformi¬ 
ties  in  the  beam  region,  and  not  from  the  presence  of  electrooptic  scattering 
centers.  The  absence  of  such  centers  is  consistent  with  results  obtained  in 
photorefractive  and  in  other  electric-field  experiments  described  in  this 
section. 
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VI.  DEPENDENCE  OF  SCATTERING  ON  OPTICAL  WAVELENGTH,  DIFFUSION 
TREATMENT,  WAVEGUIDE  VARIETY,  AND  WAVEGUIDE  SUPERSTRATE 


Owing  to  the  investigative  nature  of  this  program,  a  number  of  ave¬ 
nues  of  experimentation  were  briefly  explored  and  later  terminated  or  assigned 
to  lower  priority  status.  Although  these  avenues  were  not  followed  to  comple¬ 
tion,  they  provided  a  number  of  interesting  results  that  deserve  to  be  reported 
for  their  potential  impact  on  future  research  in  the  area  of  optical  waveguide 
scattering  reduction.  This  is  the  task  of  the  present  section.  The  work  to 
be  described  includes  interesting  experiments  to  evaluate  the  dependence  of 
scattering  on  wavelength,  diffusion  heat  treatment,  and  waveguide  superstrate. 
We  include  a  discussion  of  the  nature  of  waveguides  formed  in  LiNbO^  using  a 
low  temperature  process  employing  Ag-ion  exchange,  and  a  discussion  of  obser¬ 
vations  made  using  Nb2°5~®^ass  waveguides  provided  by  the  Air  Force. 


DEPENDENCE  OF  IN-PLANE 
SCATTERING  ON  WAVELENGTH 


The  wavelength  dependence  of  in-plane  scattering  during  this  program 
was  measured  at  488  nm,  515  nm  and  633  nm.  In  addition  semi-quantitative  ob¬ 
servations  were  made  at  1.06  ym.  Most  detailed  measurements  were  made  using 

O 

a  waveguide  formed  by  diffusing  a  720  A  Ti  film  in  flowing  O2  at  1050°C  for 
6  h.  Subsequently,  the  waveguide  performance  was  improved  (but  not  optimized) 
by  4  h  polishing.  In  all  likelihood,  a  major  source  of  scattering  in  this  wave¬ 
guide  is  Ll-Ti-0  compound  surface  contamination  that  has  not  been  totally  re¬ 
moved  by  polishing. 

Figures  20  and  21  show  the  in-plane  scattered  energy  distribution  of 
sample  88-4  at  633  nm  and  515  nm  respectively.  A  comparison  of  the  two  results 
shows  that  scattering  at  the  shorter  wavelength  is  approximately  6  dB  worse 
throughout  the  angular  range  examined.  This  indicates  a  X  ^  variation,  a 
variation  that  is  in  contrast  to  any  that  may  be  predicated  on  the  basis  of 

the  theory  presented  in  Sec.  II.  According  to  that  theory,  the  wavelength  de- 

-1  -3 

pendence  at  small  angles  is  either  X  for  surface  scattering  or  X  for  volume 
scattering.  This  would  correspond  to  a  0.9  dB  increase  in  scattering  in  go¬ 
ing  from  633  nm  to  515  nm,  if  surface  effects  dominate,  or  a  2.7  dB  increase 
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in  scattering  if  volume  effects  dominate.  The  latter  value  is  within  3.3  dB 
of  what  we  actually  observe.  Later  in  the  program  we  noted  that  the  back¬ 
ground  scattering  level  in  a  waveguide  could  change  by  about  3  dB  depending 
on  experimental  conditions,  so  perhaps  our  results  can  be  judged  consistent 

with  theory  if  volume  scattering  is  assumed.  However,  our  theory  also  states 

2 

that  the  wavelength  dependence  of  scattering  changes  by  a  factor  X  at  angles 

<(>  for  which  2irnfl/X>>l,  where  a  is  the  size  of  the  scattering  center.  For 

-3  -1 

volume  scattering  the  wavelength  dependence  changes  from  X  to  X  ,  and  the 
increase  in  scattering  in  going  from  633  to  515  nm  wavelength  should  only  be 
0.9  dB,  in  contrast  to  6  dB  observed. 

To  further  complicate  matters  we  found  that  the  in-plane  scattered 
energy  distribution  measured  at  488  nm  was  very  similar  to  that  at  515  nm.  If 
anything,  scattering  was  reduced  in  going  to  the  shorter  wavelength.  Accord¬ 
ing  to  the  measured  X  ^  dependence,  the  scattering  level  should  have  increased 
by  1.4  dB. 

These  results  are  not  supportive  of  the  theory  developed  in  Sec.  II, 
but  before  that  theory  is  modified  it  is  necessary  to  take  more  numerous  and 
more  careful  data  in  the  wavelength  dependence  of  scattering  than  we  did  in 
the  program.  As  an  example  of  the  care  required,  we  note  that  the  scattering 
associated  with  prism  coupling  may  well  be  A-dependent.  This  possibility  has 
not  been  considered  in  our  interpretation  of  the  data  of  Fig.  20  and  21.  More¬ 
over  these  data  could  have  been  influenced  by  an  experimental  problem  that  re¬ 
sulted  in  spurious  scattering  distributions  for  a  period  of  time  during  this 
program. 

DEPENDENCE  OF  IN-PLANE  SCATTERING  ON 
DIFFUSION  HEAT-TREATMENT 

In  Fig.  7  we  presented  the  results  of  a  study  that  showed  the  influ¬ 
ence  of  Ti-film  thickness  on  the  scattering  characteristics  of  TirLiNbO^  wave¬ 
guides.  These  results  motivate  the  question  as  to  whether,  for  constant  Ti- 
film  thickness,  waveguide  quality  depends  on  diffusion  time,  temperature,  or 
crystal  orientation.  Our  results  based  on  a  few  waveguide  samples  indicate 
that  this  is  not  the  case: 

O 

Item  1.  Sample  88-5  (see  Figs.  14  and  15),  made  by  diffusing  280  A 
of  Ti  for  3  h  at  950°C,  and  subsequently  polished  for  1  h,  was  given  an  addi¬ 
tional  3  h  heat  treatment.  There  was  no  significant  change  in  the  in-plane 
scattered-energy  distribution  as  a  result. 
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Item  2.  y-cut  sample  109  (see  Figs.  5  and  6)  exhibited  scattering 
characteristics  similar  to  those  of  x-cut  sample  88-5.  Both  waveguides  were 
made  under  similar  diffusion  conditions. 

O 

Item  3.  Sample  88-3,  made  by  diffusing  720  A  of  Ti  for  8.5  h  at 
950°C,  exhibited  scattering  characteristics  similar  to  those  of  sample  88-4 
(Figs.  20  and  21)  made  by  diffusing  720  A  of  Ti  for  6  h  at  1050°C.  For  the 
comparison,  neither  sample  was  polished.  Both  samples  were  x-cut  and  performed 
similarly  to  y-cut  sample  70,  which  was  comparable  to  88-3  and  4  except  for 
a  shorter  3  h  diffusion  time.  (Various  results  for  sample  70  are  displayed  in 
Figs.  7,  10,  11  and  12). 

These  items  indicate  that  waveguide  scattering  performance  is  rela¬ 
tively  insensitive  to  the  time  and  temperature  of  the  diffusion  treatment  and 
to  the  crystal  orientation  (x-cut  or  y-cut)  of  the  LiNbO^  surface.  Accordingly, 
we  did  not  devote  much  effort  to  categorizing  the  small  dependences  on  these 
conditions  that  may  exist. 

However,  at  one  point  in  the  program  we  were  concerned  about  the 
possibility  that  the  high  temperatures  of  diffusion  could  introduce  scattering 
centers  that  might  be  avoided  through  the  use  of  a  low-temperature  diffusion 
process.  In  an  attempt  to  check  this  out  we  fabricated  and  tested  an  AgrLiNbO^ 
waveguide.  Our  first  attempt  to  make  the  waveguide  involved  the  thermal  diffu¬ 
sion  of  Ag-films  at  temperatures  up  to  450°C  for  times  up  to  14  h.  In  x-cut 
sample  105,  the  film  appeared  to  be  partially  diffused.  Residue  was  nonuni- 
formly  distributed  across  the  substrate  surface.  No  waveguiding  was  observed, 
even  after  the  residue  was  lightly  polished. 

We  finally  were  successful  in  making  the  desired  waveguide  using 
Ag-Li  ion-exchange  in  a  bath  of  molten  AgNO^.^^  The  exchange  took  place  at 
a  temperature  of  360°C  for  a  time  of  6.5  h.  The  waveguide  exhibited  a  single 
TEo  mode,  and  that  was  for  the  case  of  propagation  perpendicular  to  the  optic 
axis. 

The  surface  of  the  waveguide  was  strongly  degraded  by  the  fabrica¬ 
tion  process.  It  was  improved  but  not  restored  by  a  three  hour  polishing 
treatment.  A  500X  Nomarski  micrograph  of  a  typical  portion  of  the  polished 
surface  is  shown  in  Fig.  22. 

The  waveguide  appeared  to  be  quite  lossy;  however  it  was  difficult 
to  evaluate  the  losses  visually  because  of  poor  input-  and  output-coupling 
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Fig.  22.  Nomarski  micrograph  of  the  surface  of 


an  Ag:LiNb03  waveguide  formed  by  ion- 
exchange  ( 500X) . 
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efficiency.  The  sample  had  a  slightly  brownish  cast,  and  it  is  possible  that 
much  of  the  loss  was  from  absorption  rather  than  scattering.  This  was  further 
indicated  by  the  surprising  result  that  the  scattered  energy  distribution  of 
the  AgtLiNbO^  waveguide,  shown  in  Fig.  23  was  only  slightly  inferior  to  those 
of  good  Ti:LiNbO^  waveguides.  This  result  shows  an  insensitivity  of  in-plane 
scattering  to  surface  condition  that  was,  to  us,  entirely  unexpected  on  the 
basis  of  beliefs  we  held  at  the  start  of  the  program. 

DEPENDENCE  OF  IN-PLANE  SCATTERING 
ON  WAVEGUIDE  SUPERSTRATE 


In  the  theory  of  optical-waveguide  scattering  presented  in  Sec.  II, 
a  number  of  formulas  were  derived  which  showed  a  dependence  on  the  supers trate 
refractive  index.  These  formulas  do  not  include  those  of  Eq.  (28)  for  the 
in-plane  differential  scattering  cross  section.  Consequently,  we  expect  the 
in-plane  scattered  energy  distribution  to  remain  unchanged  when  the  waveguide 
superstrate  is  changed,  as  may  be  done  by  adding  a  liquid  to  the  waveguide 
surface.  We  verified  experimentally  that  this  was  the  case,  using  a  liquid 
having  nQ  =  1.544,  obtained  from  a  set  of  index-matching  liquids. 

Among  the  formulas  derived  in  Sec.  II  which  do  depend  on  the  super¬ 
strate  index  are  Eq.  (37)  for  the  out-of-plane  differential  scattering  cross 
section  and  Eq.  (40)  for  the  waveguide  attenuation.  Assuming  that  either  sur¬ 
face  or  volume  scattering  makes  the  dominant  contribution  to  total  out-of-plane 
scattering,  we  may  write 


a 


(49) 


where  a  is  the  attenuation  coefficient,  K  is  a  constant  that  depends  on  scat¬ 
tering  and  waveguide  parameters,  n^  is  the  waveguide  index,  nQ  is  the  super¬ 
strate  index,  and  m  =  +0  for  surface  scattering  or  m  =  -3  for  volume  scatter¬ 
ing. 

"■oL 

The  quantity  P  =  PQe  is  the  power  transmitted  through  a  length  L 
of  waveguide  when  Pq  is  the  initial  power.  The  change  in  transmitted  power 
that  results  when  a  is  changed  is  an  experimentally  accessible  quantity.  In 
particular , 
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AP/P  -  e(aD  "  “w)L  -1 


(50) 


is  the  fractional  change  in  transmitted  power  that  occurs  when  a  liquid  is 
added  to  the  waveguide  surface,  where  is  the  attenuation  coefficient  for 
a  wet  surface,  and  aD, is  the  attenuation  for  a  dry  surface  (n^  =  1).  For 
low  loss  waveguides, 

AP/P  =  aDL(l  -  dy/op)  (51) 

where,  using  Eq.  49, 

ay/ap  =  [<nj  -  n*)/(n*  -  1) ]m  ,  (52) 

For  the  values  n,  =  2.2  and  n  =  1.544,  we  find 
1  o 

(4P/P)SUrf  "  0  (53> 

(SP/P)vol  *'2-82  l‘BL 

for  surface  and  volume  scattering,  respectively. 

We  have  measured  AP/P  for  Ti:LiNbo^  waveguide  sample  88-3  and  for 
AgiLiNbO^  waveguide  sample  105.  Both  waveguides  exhibited  approximately  30  dB 
scatter  reduction  at  1°,  but  the  latter  had  a  much  inferior  surface  quality, 
as  shown  in  Fig.  22.  For  the  Ti:LiNbO^  waveguide  we  find  AP/P  *  +  0.0105  + 
0.0078,  while  for  the  Ag:LiNb03  waveguide  our  measurement  yielded 
AP/P  *  0.1148  +  0.0088. 

The  fact  that  positive  values  are  obtained  suggests  that  surface 
roughness  is  the  dominant  out-of-plane  scattering  mechanism.  However,  the 
theory  of  Section  II  is  not  adequate  to  predict  these  values.  This  may  be 
owing  to  several  simplifying  approximations  used  in  the  development  of  that 
theory.  A  useful  set  of  measurements  would  be  one  in  which  AP/P  is  measured 
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for  a  range  of  superstate  indices  nQ,  so  that  the  functional  dependence 
of  a  on  nQ  could  be  accurately  tested.  Possibly,  this  technique  could  be 
helpful  in  assessing  the  relative  amounts  of  surface  roughness  and  volume 
scattering. 


IN-PLANE  SCATTERING  CHARACTERISTICS  OF  AN 
Nb205-GLASS  WAVEGUIDE 

Alternatives  to  Ti:LiNbO_  waveguides  that  are  candidates  for  use  in 
the  spectrum-analyzer  application^  are  waveguides  formed  by  sputtering  glass 
or  metal  oxides  on  thermally  oxidized  Si  substrates.  We  report  briefly  here 
on  observations  we  made  on  an  Nb20^-microscope-slide  waveguide  that  was  pro¬ 
vided  us  by  Dr.  Douglas  Wille  of  the  Air  Force  Avionics  Lab. 

The  in-plane  scattered-energy  distribution  for  this  waveguide  is 
shown  in  Fig.  24.  The  asymmetry  in  the  angular  distribution  is  believed  to 
result  from  a  geometric  asymmetry  in  our  detection  system  that  we  were  experi¬ 
encing  at  the  time  this  graph  was  generated.  Consequently,  repeat  data  should 
be  acquired  before  final  evaluation  of  this  waveguide  relative  to  LiNbO^  wave¬ 
guides  tested  in  the  program.  Our  tentative  conclusion  is  that  the  two  types 
are  comparable  with  regard  to  the  initial  rate  of  fall-off  of  scattered  in¬ 
tensity  with  angle.  However,  the  total  angular  span  of  the  visible  in-line 
was  much  greater  in  the  Nb20^-glass  waveguide  than  in  a  typical  Ti:LiNbO^ 
waveguide.  The  longer  m-line  would  be  expected  to  result  from  the  presence 
of  a  background  distribution  of  particles  small  in  size  compared  with  a  wave¬ 
length.  The  apparent  absence  of  such  a  distribution  in  LiNbO^  waveguides  is 
an  attractive  feature  of  this  material. 
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